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Abstract. One of the most interesting sources of gravitational waves (GWs) for LISA 
is the inspiral of compact objects on to a massive black hole (MBH), commonly referred 
to as an "extreme- mass ratio inspiral" (EMRI). The small object, typically a stellar 
black hole (bh), emits significant amounts of GW along each orbit in the detector 
bandwidth. The slowly, adiabatic inspiral of these sources will allow us to map space- 
time around MBHs in detail, as well as to test our current conception of gravitation 
in the strong regime. The event rate of this kind of source has been addressed many 
times in the literature and the numbers reported fluctuate by orders of magnitude. 
On the other hand, recent observations of the Galactic center revealed a dearth of 
giant stars inside the inner parsec relative to the numbers theoretically expected for 
a fully relaxed stellar cusp. The possibility of unrelaxed nuclei (or, equivalently, with 
no or only a very shallow cusp, or core) adds substantial uncertainty to the estimates. 
Having this timely question in mind, we run a significant number of direct-summation 
A— body simulations with up to half a million particles to calibrate a much faster orbit- 
averaged Fokker-Planck code. We show that, under quite generic initial conditions, the 
time required for the growth of a relaxed, mass segregated stellar cusp is shorter than 
a Hubble time for MBHs with M. < 5 x 10 6 M Q (i.e. nuclei in the range of LISA). 
We then investigate the regime of strong mass segregation (SMS) for models with two 
different stellar mass components. Given the most recent stellar mass normalization for 
the inner parsec of the Galactic center, SMS has the significant impact of boosting the 
EMRI rates by a factor of ~ 10 in comparison to what would result from a 7/4— Bahcall 
& Wolf cusp resulting in ~ 250 events per Gyr per Milky Way type galaxy. Such 
intrinsic rate should translate roughly into ~ 10 2 — 7 x 10 2 sbh's EMRIs detected by 
LISA over a mission lifetime of 2 or 5 years, respectively), depending on the detailed 
assumptions regarding LISA detection capabilities. 
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1. Introduction 



Nowadays it is well-established that a massive dark object, very possibly a massive black 
hole (MBH) with a mass of about 4 x 10 6 M^, is lurking in the centr e of the Milky Way 
jEisenhauer et alihoosi bhez et al.lhoO.j boosi loillessen et aDhood ). While there is an 
emerging consensus about the origin and growth of supermassiye black holes (SMBH , 



with masses ab out or larger than 1O 8 M ) (IFerrarese fc Fordll2005l ; iDiMatteo et al.ll2008 



Volonterill20inL MBHs with smaller masses such as the one in the Galactic centre remain 
a (relatively) understudied enigma. One of the keys to understanding the growth and 
evolution of MBHs in this lower mass range resides in the dynamics of stars in their 
vicinity. This is the case mainly because relaxation times there are low enough that the 
surrounding stellar systems should have had enough time — through two-body relaxation 
alone — to evolve towards a steady-state which is independent of the particular initial 
conditions at the time of formation. The Galactic center is thought to fulfill such 
condition. It is the universality of such relaxed stellar nuclei that gives us a crucial 
predictive power on the expected properties of the MBH environment, on the stellar 
candidates for close interaction with the central MBH and on the resulting gravitational 
wave (GW) signatures. If, on the contrary, non-relaxed systems were generic, then one 
would need to resort to case- by- case modelling of each galactic nucleus. 

The ideal probe for these innermost regions of galaxies is the GW radiation that is 
emitted by stellar bhs and other compact objects that come very close to the MBH. One 
of the main channels for interaction between stars and a central MBH is the adiabatic, 
slow inspiral of compact remnants (CR) into the MBH due to the emission of GWs — an 
EMRI. During such an event, the small body effectively acts as a probe of spacetime 
close to the MBH as its orbit slowly shrinks due to the energy and angular momentum 
lost in the form of GW radiation. In case of 10 5 — 1O 6 M MBHs, after some ~ 10 4 — 10 5 
orbits in the LISA band (/ or b > 1CT 4 Hz and a periapsis a < few x Rs c hw, since we 
only consider sources which are completely embedded in the band, and not bursting 
sources), the small body eventually merges with the MBH. The information contained 
in the waves will allow us to determine t he parameters of these binary system with a n 
unprecedented accuracy (see for instance iBarack and Cutlerl l2004i ; iBabak et al.ll2010l ). 
corroborate the existence of MBHs and maybe even provide the first direct detection of 
an intermediate MBH (in case the primary is ~ 1O 3_4 M ). 

LISA will thus scrutinize exactly the mass range about which electromagnetic 
observational information is currently lacking. In its most general form, the EMRI 
problem — the astrophysical modelling of event rates and parameters for EMRIs — 
spans many orders of magnitude. From the bulge regions at fewxlO pc, where the 
dynamics is essentially collisionless -but from where single stellar bhs and binaries with 
CRs originate; down to the parsec scale of the nucleus itself which evolves secularly 
over (local) relaxation timescales; and then further down to milliparsec scales where 
relativistic effects start to dominate the evolution. But, however, once a steady state 
configuration establishes itself in the central parsec region, the EMRI rates are rather 
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expected to depend strongly on the (universal) density distribution of CRs within (in 
order of magnitude) O(O.Olpc) from the hole. This is indeed the region from whic h 
these inspiralling sources are expected to originate ( iHopman and Alexander! 120051 ). 
The dynamics in this tiny volume has been rather unexplored until the relevance of 
EMRIs and sub-parsec observations of the Galactic center have raised its interest. 
Since then, ma ny authors have devoted a number of works to the analysis of this 



peculiar regime (pigurdsson and Rees 



1997; 



Freitag 



2003 



Alexander and Hopman 



2003 



Hopman and Alexander! 120061 ) . 



We discuss in this work the stellar distribution of dense stellar systems 
around MBHs in the LISA mass range. Realistic modeling of mass segregation — 
which is the natural outcome for any realistic stellar population — will strongly 
impact the expected E MRI rates, since it favors the accumulation of heavier 
objects towards the center ([Hopman and Alexanderll2006l : lAlexander and Hopmanll2009l ; 



Preto and Amaro-Seoand l2010al). In S e ction 2 , we begin by summarizing the results 



obtained by iPreto and Amaro-Seoand ( l2010al ) that show how to calibrate the FP 
calculations with direct A^-body simulations^]; then, still in the same section, we present 
new results concerning the robustness of A^-body realisations of stellar cusp growth with 
respect to the total particle number N. In section 3, we present new results on the growth 
of stellar cusps from a variety of initial conditions resulting from carving a cavity in the 
star's phase space distribution function. This is motivated from a variety of astrophysical 
mechanisms that may lead to cusp destruction; and these mechanisms are critically 
assessed in the end of the section. We show that, under very generic circumstances, 
the time required for the growth of a cusp is shorter than a Hubble time. Therefore, 
quasi-steady, mass segregated, st ellar cusps are expe cted to be common a round MBHs 
in the LISA mass range. But see iMerrittl ( 120101 ) and lMadigan et al.l (j2010[ ) for different 
perspective s. EMRI detection rates for LISA are expected to peak for M, ~ 10 5 — 1O 6 M 
(lGairll2009i ) leading us to conclude that at least a sizeable fraction of these events should 
originate from strongly segregated cusps. Finally, in Section 4 we present new estimates 
on the expected EMRI rates in mass segregated nuclei and conclude that our realistic 
modeling of mass segregation has a significant impact on these rates. 



2. Mass segregation 



The distribution of stars around a massive black hole is a classical problem in stellar 



dynam ics ( Bahcall and Wolf Il976t iLightman and Shapiro! Il977l ). iBahcall and Woll 
Jl976h have shown, through a kinetic treatment that, within the radius of gravitational 
influence of the hole r^, in case all stars are of the same mass, this quasi-steady 
distribution takes the form of power laws, p(r) ~ r~ 7 , in physical space with p(r) 
the stellar density at a radius r and f(E) ~ E p in energy space (with E the 

| Direct iV-body simulations compute the gravitational accelerations between particles using the 
exact Newton's law, withou t introducing any approximations to compute the gravitational potential 
( Binnev and Tremaine 20081) . 



Event rate of EMRIs and cusp re-growth 



energy and 7 = 7/4 and p = 7 — 3/2 = l/4)[§. This is the so-call ed zero-flow 
soluti on for which the net flux of st ars in energy space is precisely zero. IPreto et al 



( 12004J ) and iBaumgardt et al.l (120041 ) were the first to demonstrate the robustness of 
the corresponding direct-summation iV-body realizations, and have therefore validated 
the assumptions inherent to the Fokker-Planck (FP) approximation — namely, that 
scattering is dominated by uncorrelated, 2-body encounters and, in particular, dense 
stellar cuspsjjj] populated with stars of the same mass are robust against ejection of stars 
from the cusp. The latter result is not trivial as for a BW 7 = 7/4 cusp stellar densities 
are extremely high at the center and the fraction of stars with speeds close to the escape 
velocity fr om the cusp is quite high at all radii r < r^, with r^ the influence radius of 
the MBH JPretoll2010h . 

Single mass models are very poor approximations of real stellar populations. To 
first order of approximation, an evolved stellar population can be represented by two 
(well-separated) mass scales: one in the range 0(1M Q ) corresponding to low mass main- 
sequence stars, white dwarfs (WDs) and neutron stars (NSs); another with (9(1OM ) 
representing stellar bhs. Therefore, for simplicity, here we restrict our discussion to 
models with two mass components and leave the more general case to another work in 
preparation (IPreto and Amaro-Seoandl2010bl ). 

When stars of two different masses are present, there is mass segregation which 
is a process by wh ich the heavy stars accumulate near the center while the lighter 



ones float outward (jSpitzer 



1987 



Khalisi et al.l 120071). Accordingly, stars w ith different 



mass get distributed with different density profiles. iBahcall and Woffl (119771 ). henceforth 
BW77, have argued heuristically that a scaling relation p^ = mi/mj x pj (where the 
sub indices i, j refer to the light or heavy components) establishes itself and depends 
only on the mass ratio. Here, as in the single-mass case, the crucial assumption is 
that all components are abundant enough that they undergo enough scattering among 
themselves and with the other components as to stabilize into an approximate zero-flow 
solution. Obviously, t his cannot happen independently of the number fraction of th e 
different stellar masses (lAlexander and Hopmanll2009l : IPreto and Amaro-Seoane!l2010al ) . 
In the realistic situation where the num ber fraction of heavy obje c ts (in our case, stellar 
bhs) is small, a new solution coined by lAlexander and Hopmanl (120091 ) as strong mass 



segregation (SMS) obtains with density of heavy objects scaling as Ph(j) ~ r 
a > 2. The solution has two branches and can be parametrized by the parameter 



where 



A 



D (l) +D 

HH + u 



(2) 
HH 



N H m 2 H 



D (1) + D 

LH + u 



(2) 
LH 



Ni,m L 3 + m/f/mi' 



(1) 



where Nl and Nh are the total number of light and heavy stars, mi and m# are the 
corresponding individual masses. A provides a measure of the importance of the heavy 



§ We note that 12 years before the work of BW, iGurevich ( 1964 ) derived a similar solution for how 
electrons distribute around a positively charged Coulomb center, which is the equivalent of the MBH 
in our case. 

|| In this work, a nucleus is said to be a core if 7 < 1; it is said to be a cusp if 7 > 1. 
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star's self-coupling relative to the light-heavy coupling (in terms of the 1 st and 2 nd order 
diffusion coefficients); and it depends essentially on the mass and number ratios, which is 
one parameter more than proposed by BW77. The weak branch, for A > f corresponds 
to the scaling relations found by BW77; while the strong branch, for A < 1, generalizes 
the BW77 solutiorfl. Stellar populations with continuous star formation and an initial 
mass function (IMF) given by dN/dM oc M~ a will be characterized by A < 1 if a > 1.8 
and A < 1 otherwise; and, in p articular, Salpeter and Kroup a's IMF generate evolved 
stellar populations with A < 1 ([Alexander and Hopmanll2009l ). 

Figure 1. Evolution of density profiles. Mass density profiles, p^(r) (left panels) 
and pH{f) (right panels) at the end of the integrations, after « 0.2T r i x (rh). Red curves 
are from FP calculations, green 
curves are from NB simulations. The agreement 
between both methods is quite good. The mass 
ratio between heavy and light stars is R — 
10, representing the expected typical mass ratio 
between light stars (MS stars, WDs and NSs) and 
stellar bhs, as explained in the text; the number 
fraction of heavy stars /// = 2.5 x 10 -3 (top and 
middle panels) and fn = 0.429 (lower panels), 
corresponding to the strong and weak segregation 
regimes respectively. The initial condition is a 
Dehnen profile with ce ntral slope 7 = 1 for th e 




top an d bottom panels (jPreto and Amaro-Seoane 
2010aj), 7 = 1/2 in the middle panel; a central 
MBH with 5% of the total mass of the cluster and 
1% likewise. The particle number is N = 124, 000 
(top and bottom) and N = 512,000 (middle). The A " 3 V ° 1 " "r'' ° 1 

asymptotic slope 7# decreases from > 2 to ~ 7/4 when moving from the strong to the 
weak branch of the solution. The asymptotic slope jjj ~ 3/2 throughout, or just slightly 
below this value. The arrows point to radii r^ and O.lr^. 

There is a straightforward physical interpretation for the strong branch of mass 
segregation. In the limit where heavy stars are very scarce, they barely interact with 
each other and instead sink to the center due to dynamical friction against the sea of 
light stars. Therefore, a quasi-steady state develops in which the heavy star's current is 
not nearly zero and thus the BW77 solution does not hold exactly anymore. Indeed, in 
the limit where the number fraction /# of heavy stars is vanishingly small, as is the case 
of nuclei with realistic IMFs, the stellar potential is dominated by the light component. 
In this case, the light stars should evolve as if in isolation and develop a 7^ ~ 7/4 density 
cusp. The scarce heavy stars sink to the center due to dynamical friction against the 

% The choice of the names is based upon the resulting slopes in the density profiles, which are steeper 
{stronger) or shallower {weaker) 
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backg round of light stars, and will not exert any significant back-reaction on them (IPreto 



20101 ). 

Figure 1 displays the FP and NB evolutions of the spatial density pz,(r) and pj/(r) 
for models with two mass components corresponding to different initial profiles, MBH 
masses and total particle number N. The starting models are either 7 = 1 or 7=1/2 
Dehnen profiles for both components with a MBH of 1% or 5% of the total mass of the 
cluster. Dehnen density profiles are defined by p(r) = (3— ^)Mtot / { r b J r r Y~ 1 1 have 
total stellar mass Mtot, an inner (outer) logarithmic slope 7 (4 — 7) and a break radius 
rt> = 1 (which is larger than r% in all cases). We adopt units where G = Mtot — 1- The 
density of both components reaches a quasi-steady state wi thin ~ 0.2T r ^(rh), where 
T r h(rh ) relaxation time measured at the influence radius (IPreto and Amaro-Seoane 
2010ah . The top and middle panels display the strong mass segregation solution with 



■jh ~ 2.1 as expected for A = 0.08 (fn = 2.5 x 10~ 3 ); while, in the bottom panel, 
A = 13.2 (fx = 0.429) displays the weak solution for which 7^ ~ 7/4. The former 
value was chosen to be close to the number fraction of stellar bhs to be expected from 
a standard Salpeter or Kroupa's IMF; the latter value is chosen to be representative of 
the regime of weak segregation studied by BW77. One can see from Figure 1 that in 
the case of weak segregation Ph > Pl everywhere due to the extremely high number of 
heavy objects; in contrast, in the SMS regime pu > Pl only for r < 0.01r h (and the light 
objects dominate in number almost everywhere). In all cases the asymptotic slopes are 
valid within ~ O.lr^ and are totally predictable onc e A is known. These results ag ree 
pretty well with the predictions for the SMS regime ([Alexander and Hopmanll2009l ). 

The particle number in our direct-summation A^— body simulations sample ranges 
from A" = 124, 000 to A^ = 512, 000; our results do not show evidence of any dependence 
on total N, nor on the initial value of 7, once the results are re-scaled appropriately 
(i.e. measured in terms of the relaxation time). The agreement between NB and FP 
methods is quite good in all cases. 



3. Cusp Re-growth 

3.1. Current observations: A missing cusp 

We have seen that theory predicts a steady state cusp that reaches extremely high 
densities in the center near the MBH. Furthermore, given a normalization at, say, 
and a knowledge of the stellar mass function (and thus of A), the density profile inside 
rh becomes completely determined. But observations are much more complicated to 
interpret. First, one must realize that there are very few galaxies for which the influence 
radius can be resolved. In fact, except for the nearest galaxies, covers an angular 
region in the sky which is too small to be resolved even with the HST. The HST has an 
angular resolution of ~ 0.1". In the case of galaxies in the Virgo cluster, at ~ 17 Mpc of 
distance, it can only resolve regions with linear sizes > 8.25 pc. Therefore, HST would 
not be able to resolve SgrA*'s radius of influence if it were at the distance of Virgo. 
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Since oc MJ , it can only start to resolve the influence radius of Virgo's MBHs that 
have masses M. > 4 x 10 7 M©. Second, even in the few cases for which r^ can be 
resolved to some extent, it still is necessary to assess whether the observed stars (only 
those that are bright enough to be detected) really trace the underlying (dynamically 
dominant) invisible population. Third, given the fact that, as we have seen, stars tend 
to segregate by mass, there is an extra uncertainty related to the unknown stellar mass 
function. Moreover, there are indications that star formation events are common in 
galactic nuclei and furthermore that the resulting IMF in these sub-parsec regions may 
be substantially dif ferent from that of the field stars and biased towards heavy masses 
( IBartko et al.ll2010l ). Finally, it is necessary to deproject the observations and, in the 
(inevitable) absence of complete knowledge of phase space coordinates, one must rely 
on kinematic assumptions regarding the (an-)isotropy of stellar velocities and on the 
three dimensional shape of the stellar system. 

Figure 2. Time for cusp re-growth. Single-mass relaxation time at r^ for single- 
mass cored models as a function of MBH mass. The shaded area covers 
[0.1T rlx , 0.2T rlx ] — the time for cusp re- 
growth if there is no hole in the initial 
DF. The three dashed lines above the 
shaded region represent the average 
time needed for the cusp re-growth 1 
in case one imposes an initial cavity 3 
with size i? cav = 0.5, 1 or 2 pc. The 
horizontal dashed curve represents 
13 Gyr. It can be seen that the 
time needed to re-grow a cusp around 
MBH with masses M. < 5 x 1O 6 M 
is below a Hubble time so long as the initial cavity is smaller than < 2 pc. 

Nevertheless, it has come as a surprise that very recent spectroscopic observations 
of the Galactic center revealed a core (or even a dip) in the surface distribution of the 
old stellar population (e ssentially red giants) which should ha ve had time to relax into 
a cuspy density profile (IBuchholz et al.l 120091 ; iDo et al.l 120091 1 . The caveat is that the 
detected stars are still a small fraction, of about 5%, of the stellar population as a whole 
and therefore do not exclude the presence of an extended dark cluster (presumably 
made of stellar bhs and other CRs) — which would indeed agree with our theoretical 
expectations. 




I°9ia M BH (M Sm ) 



3.2. Carving a hole in the stellar distribution 

To assess the likelihood that the Galactic center is indeed unrelaxed, it is natural to 
ask: how long does it take to re-grow a cusp if, at some point, it has been destroyed? A 
complete answer will, of course, depend on the extent to which the cusp was destroyed, 
i.e. how much mass was expelled from the original cusp and over which radial range. 
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At this level, it does not matter very much which mechanism led to the destruction of 
the cusp. We discuss briefly possible scenarios for cusp destruction at the end of the 
section. 

In order to investigate this question, we have concocted a set of initial conditions 
purported to mimic the outcome of a destroying cusp event — such as the carving of a 
cavity in phase space through the ejection of stars by, say, an infalling IMBH or, following 
a major merger, by a MBH. We model the outcome of such an event by imposing that all 
stars with binding energies larger than some Eq or almost equivalently, with semimajor 
axis smaller than GM 9 /2E Q , are not present in the initial DF. In fact, inside r^ the 
MBH dominates the gravitational potential and E ~ GM,/2a. We thus set up an initial 
Dehnen model with f(E) = for E > E — in other words, there is an initial cavity in 
the phase space DF, but not in physical space as the stars with lower energy still entail 
p{r) ~ r -1 / 2 at the center, although with a smaller amplitude than the original model. 
The values of E were chosen to represent cavities of size R c = 0.5, 1, 2 pc resulting from 
the inspiral of an IMBH/MBH. Note that these models are, by construction, isotropic 
in the velocity distribution. El Our fiducial model is a Milky Way type nucleus with 
M. = 4 x 1O 6 M , some 1O 6 M in total stellar mass inside 1 pc distributed according to 
an initial central density slope 7 = 1/2, two components with masses = 1M and 
ran — 1OM , and 0.1% of stellar bhs by number. When the stellar distribution has no 
phase space cavity, this translates into having stars down to roughly 10~ 5 pc. Having 
validated the FP models with detailed iV-body simulations, we now proceed in the rest 
of the paper to describe results obtained with the (much faster) FP approach. 

Figure 2 shows the times for cusp re-growth computed with FP for different galactic 
nuclei models. The time for cusp re-growth is defined as the time it takes for a given 
initial density profile (p(r) in space or f(E) in phase space, with or without an initial 
cavity) to reach its asymptotic slope, which depends on A, down to r ~ O.Olr^. This 
i s indeed the scale w hich is resolved by recent observations of the Galactic center 
(jSchodel et al. 1 120091 ) . The shaded region represents the time of cusp re-growth for 
a range of R and fn (all in the SMS regime, A < 1) for the case where f(E) extends to 
high E without any cut. It can be seen that, for M. < 1O 7 M , cusps grow in less than 
a Hubble time; in the particular case of the Milky Way nucleus with M. ~ 4 x 1O 6 M , 
it takes no longer than ~ 4.8 Gyr to fully re-grow a steady-state, mass segregated, 
stellar cusp and only ~ 2.4 Gyr to have it grown down to 0.017V If, instead, an initial 
cavity is imposed at the center with size _R cav = 0.5,lor2pcin case of the Milky Way 
(or i? cav = 0.2r h ,0Ar h or 0.8r^ in case of a generic nucleus), times for re-growth are 
represented by the dashed curves above the shaded region. In this case, times for cusp 
re-growth increase; in the Milky Way case, it becomes ~ 4.8, 7.2 or 12 Gyr, respectively. 
Note that, in the mass range 10 5 — 1O 6 M , the times for cusp re-growth are definitely 
much shorter than a Hubble time — even if a fairly large cavity of size comparablee to 

+ We assume that the timescale for isotropization of velocities is much shorter than that associated 
with the cusp re-growth; in any event this should not affect our estimates by more than 10% or 20% 
maximum. 
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is hypothesized. The full curve represents the relaxation times computed at the radius 
of influence r^, while the dashed curves represent the actual times for cusp re-growth as 
measured from FP calculations (IPreto and Amaro-Seoanell2010al ). 

It is difficult to devise plausible mechanisms for the formation of such large cores in 
the stellar distribution. For instance, the inspiral of an IMBH of mass M. ~ 1O 3 ~ 4 M 
that forms an unequal-mass binary with the MBH and ejects stars through three body 
encounters would tend to progressively wipe out the stellar cusp. However , the core 
radius carved by such an event is r c ~ 0.02 — 0.04 pc (IBaumgardt et al.ll2006l ) and thus 
a steady inflow of such IMBHs (one every 10 7 years for a Hubble time) would be required 
in order to c arve a large core 50 or 100 l arger. Such large inflow of IMBHs have been 
proposed by iPortegies Zwart et al.l (120061 ). This does not seem very likely anymore in 
light of the fact that such IMBHs were hypothesized to be formed by runaway mergers of 
stars in the center of globular clusters. However, at solar metallicities, such mechanism 
seems very inneficient. Mass loss due to very strong winds severely limits the growth 
of the stellar object being formed and the likely end result of a runaway merger is a 
~ IOOMq Wolf-Rayet star. At lower metallicities, mass loss is lower an d the remnant 
can be more massive ~ 26OM , but in any case it will not form an IMBH (IGlebbek et al. 



20091 ). In sum, it looks very unlikely that sufficient IMBHs can be formed in order to 
generate such steady inflow to the Galactic center. Another possibility would be that 
SgrA* is a binary MBH, but this would most likely imply that there has been a, more or 
less recent, major merger involving the Milky Way. This would contradict the apparent 
pure-disk nature of the Galaxy, as theoretical interpretations of stellar kinematic data 
of the Galactic Bulge seem to favor that the Bu l ge is p art of the disk and not a separate 
component resulting from a merger IShen et al.l (]2010l ) y — aside from the fact that there 
are strong constraints from the SgrA* proper motion (IReid and Brunthalerl 120041 ) . 

Stars in a Keplerian potential, GM,/r, do not precess because of the 1 : 1 resonance 
between their radial and azimuthal frequencies. Resonant relaxation (RR) results from 
the coherent torques that such stars exert on each other leading to a fast evolution 
of their orbital angular momenta over timescales < T pr , where T pr is the precessio n 



timescale due t o depa rtures from an exact Kepler potential (IRauch and Tremaindll996l ). 
Madigan et al.l (120101 ) suggest that RR, by increasing the rate of angular momentum 
diffusion in the near-Keplerian gravitational potential around the MBH, may boost the 
tidal disruption rate of stars and could thus create a near-cavity (out to ~ 0.1 pc) 
in the stellar distribution. It is certain that RR operates to some extent in the inner 
parsec, but we doubt it can completetly explain the dearth of red giant stars there or, 
more generally, the full destruction of a cusp — including CRs such as stellar bhs3 First, 
their final density distribution does not show a cavity, nor a shallow cusp profile, for 
r > 0.1 pc; instead, they get final slopes 7 ~ 1.5 for r > 0.1 pc. Th is is in contrast 



with observations which show a decaying density for r < 0.24 pc ( jBuchholz et al 



20091 ). Second, we believe that their Monte-Carlo calculations severely underestimate 



* Resonant relaxation is. never theless, very likely to have a significant impact on the EMRI event rates 
( Hopman and Alexander 20061) . 
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the rate at which the cusp re-grows; in fact, they obtain a timescale which is ~ 10 — 30 
times longer than that obtained with our iV-body simulations (which are free from 
any simplifyi ng dynamical assumptions), and also from ours and their own recent FP 



calculations (IHopman and Madiganll2010l ). Third, were they to use the latter rates, and 
given that the time RR takes to deplete the cusp is of the same order as the time we 
obtain for cusp re-growth, the net effect of RR on the cusp would likely be minute. 
Moreover, stellar bhs cannot be tidally disrupted, and this makes them less susceptible 
to be extracted from the cusp than 1OM stars. 

4. EMRIs rates 

4-1. Adiabatic and abrupt EMRIs: Estimation of the rates 

Given a steady state stellar bhs continue to diffuse in (E, J)-space and some of them 
eventually come into close interaction with the MBH. During a close interaction, a 
stellar bh can either be promptly scattered into the MBH, accompanied by a single or 
a few brief bursts of GWs in the LISA band — the so-called "direct-plun ges" — , though 



they are not likely detectable unless if emitted from the Galactic center (IHopman et al. 



20071 ). or scattered outwards in the cusp. In either case, it does not live enough to 
become an EMRI. Alternatively, it may undergo a very slow, adiabatic, inspiral without 
being appreciably disturbed by other stars and, in this case, it will eventually become 
an EMRI detectable by LISA. An EMRI object thus has to spend very many orbits 
without being significantly scattered by the gravitational tugs of the other stars. In 
other words, they must fullfill the following inspiralling criterion: the time Tgw it takes 
for the inspiral, due to orbital energy lost by GW emission only, must be shorter than the 
typical time Tj it takes on average to drift in angular momentum by an amount J which 
equals its orbital angular momentum. Otherwise, they will be promptly captured by the 
MBH before entering the LISA band. The inspiral criterion can be stated in terms of 
the parameter s being smaller than unity, s = T GW /Tj < 1. For T GW > Tj, it is almost 
certain that this object has either taken an almost radial orbit and fallen into the MBH 
as a direct plunge or has been scattered outwards. It turns out that this paramete r 
simply scales with orbital's semimajor axis: s oc c?l 2 ~ v ( Hopman and Alexander! 120051 ). 



which means that it is a decreasing function of a so long as y < 3/2. This is indeed 



the ca se in both regimes of mass segregation. Furthermore, IHopman and Alexander 



( 120051 ) have shown that the probability for a successful inspiral as a function of orbital 
semimajor axis (or energy) is almost a step function of semimajor axis. If a < acw? 
it is almost certain that the stellar bh will become an EMRI; it will almost certainly 
not become one in case the inequality sign is reversed (and the width of the "transition 
region" is very small). This crucial threshold quantity demarcates the orbits which are 
close enough to the MBH to sucessfully decouple from the rest of the cluster and undergo 

() In steady state, on average each star that drifts outward by an amount J will be compensated by 
another that drifts inward by the same amount. This balance only breaks down for those orbits that 
fall on to the hole, since there are obviously no stars coming out of it to keep detailed balancing. 
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the slow, adiabatic inspiral that defines an EMRI from those more weakly bound orbits 
that will be perturbed out of the EMRI tracks due to scattering with other stars. 

Therefore, in order to estimate the EMRI event rate, given a steady state f(E), 
obtained via the FP equation, one essentially counts the number of stars that populate 
the region of phase space for which the inspiralling criterion above is satisfied and divide 
it by the local relaxation time. Note that here, for simplicity, we ignore other driving 
mechanisms — in particular, we ignore resonant relaxation. Under these assumptions, 
the EMRI rate for stellar bhs is approximately given by 



EMRI 



n(E) 



E GW HJc(E)/J lc ) 



(2) 



T lU (EY 

where /. is the number fraction of SBHs in the stellar popu lation, n(E) is the total 
number of stars per unit energy(n(£') oc f(E), see IPretol ( 120101 )). J C (E) = y/GM,/2E is 
the specific angular momentum of a circular orbit of energy E, J\ c = 4 GM,/c is the loss- 
cone angular momentum and T r j x = 0.34 a 3 /[G 2 (m 9 p 9 + m^p*) In A] is the relaxation 
time. The log term in the denominator in ([2]) arises from the phase space (partial) 
depletion resulting from the presence of the loss cone. The conversion between r and E 
is, for r <C r^, (E(r)) = GM,/2r or E = GM,/2a. The critical radius acw, or energy 
E GW , f or EMRIs is approximately ag w = 0.01r fe ; and, to first order, acw is independent 
of M. (jHopman and Alexander! 120051 ) . 



4-2. The relevance of realistic models of mass segregation for the rates 

The weak regime of SMS, and corresponding BW solution, would lead to a fairly high 
intrinsic rate, per galaxy, of EMRIs. In fact, Figure 3 shows that, for a Milky Way 
nucleus, in case A > 1, the intrinsic EMRI rate is > 10 3 per Gyr. This is, however, 
unrealistic as such scenario pressuposes an unrealistically high number fraction of bhs 
(/• ~ 0.0325 for A > 1). In the more realistic case, when A ~ 0.03 the BW solution 
would entail a strong supression of the EMRI rate to — at best — a few tens of events per 
Gyr. This is where SMS solution appears to rescue us. SMS implies a higher density of 
bhs inside as compared with the 7 = 7/4 solution, and in this way — by decreasing 
the local T r \ x and increasing n(E) close to the MBH — it partially, but not completely, 
compensates for the small number fraction of bhs entailed by realistic mass functions. 

Figure 3. EMRI rate as a function of A. The number of stellar bh EMRI events 
per Gyr in a Milky Way type nucleus (M. = 4 x 1O 6 M and M*(< lpc) = 1O 6 M ) 



Event rate of EMRIs and cusp re-growth 

as a function of the parameter 
A. This is computed from a 
two-component mass segregated 
stellar cusp (7^ ~ 2.1 and 
7l ~ 1.5) with mass ratios R = 
10 and 15 obtained from FP s 
calculations. In the case of the £ 
fiducial value /. = 1(T 3 , A ps I 
0.03; in those circumstances, 
each Milky Way like nucleus will 
produce on average ~ 250 stellar 
bh EMRIs per Gyr. 
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Milky Way type nucleus, M f 




log 10 A 



In order to quantitatively evaluate the boost r SMS /r B w to the EMRI rates from 
SMS, for a given A and a fixed mass normalization at r^, one needs to estimate what 
would be the rate if the spatial and phase space densities were determined by the 7 = 7/4 
cusp for r < O.lr^. This is done as follows: we define analytically both p(r) and f(E) 
that would result from a 7 = 7/4 inside O.lr^ 

p(r) = p F p(r), r > r L 

P(r) = p F p(r L ) x (-j , r <r L , (3) 

and 



f{E) = fpp(E), 

f{E) = f FP (E L ) x (J- 



E <E L 

1/4 



E>E L , 



(4) 



Milky Way type nucleus, M BH =4x1 M 



R=10 
R=15 



where the indices FP mean that the profile is taken from the Fokker-Planck calculation. 
r L (and E L ) is a reference radius (energy) chosen according to r L ~ 0.1r h . 

Figure 4. Boost on the EMRI rates from strong segregation. 
One can see that, for values of 
A < 1, there is a significant 
boost to the EMRI rates in 
comparison to which it would 
result in the case of a 7 = 7/4 
BW cusp. In particular, for our t 
fiducial value A ~ 0.03 (/. ~ ^ 
10 ), the boost is of order of a ° 
factor 10 with respect to a 7/4- 
BW cusp with the same mass 
normalization at r = 1 pc. 




Iog 10 A 



Figure 4 shows the boost to 
the EMRI rates due to SMS relative to what would be obtained from a BW profile. 
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Going from an unrealistically high /., as adopted by BW77 (say A = 3), to a more 
realistic /. (say A = 0.03), while neglecting the existence of SMS, one supresses the 
EMRI rate by factors of ~ 100 — 150 (the former would lead to ~ few x 10 3 EMRIs per 
Gyr; the latter is reduced to ~ few tens per Gyr). However, by taking into account the 
SMS solution, for this low A = 0.03, we boost again the rates by a factor close to 10, 
thus partially compensating the reduction of EMRIs (from few tens to a few hundred 
per Gyr; in fact, there are ~ 250 per Gyr in case A = 0.03 for a Milky Way nucleus). 
We conclude that the apparently inocuous and tiny change of the logarithmic slope from 
7# = 7/4 to jg ~ 2 can have a substantial effect (a factor of ~ 10) on the expected 
EMRI rate. 

Figure 5. EMRI rates as a function of MBH mass in strongly segregated nuclei. 
The EMRI rate depends on the 

-0.19 





H-10,f Bh -10| —— ' 
R=10,f BH =10* — — 
R=15, f BH = 1 — *— 


^ -0.19 — 




-0.19 





logic M BH ( in M o) 



MBH mass, T SMS oc M, 
Shown are curves for A = 0.03 
(/. = 10- 3 and 10~ 2 ) and three 
different mass ratios between 
heavy and light stars, R = 10 
and 15. 

Figure 5 shows the depen- 
dence of the intrinsic EMRI 
rate on the mass of the central 
MBH, where the validity of the 
M. — a relation was assumed 
( iFerrarese fc Fordll2005l ). We see that the EMRI rate for stellar bhs scales as oc M" 019 , 
independently of R and Its absolute normalization depends obviously on the number 
fraction /. of sbhs, in agreement with Figures 3 and 4. 

One can ma ke a rough co nversion of the estimated intrinsic rate into LISA detection 
rates. Following iGairl ( 120091 ). who made a number of assumptions regarding the local 
density of MBHs and its spin distribution, plus on the LISA detection capabilities, 
we find that according to its equation (7), LISA will see around ~ 10 2 — 7 x 10 2 
EMRI events during a 2-year or 5-year mission, respectively. Note that these rates 
may cha nge by factors of ~ 2 — 3 as a function of corrections to the local MBH mass 
density (IGraham and Driverll2007l ); moreover, the uncertainties regarding the efficiency 
of RR and other channels may still affect the rate predictions by one or two orders 
of magnitude. A lot of work still remains to be done; nevertheless, the consequences 
regarding the SMS regime are significant and under control. 



5. Conclusions 



We have considered simplified stellar models of galactic nuclei, with only two mass 
components, which harbor MBHs that fall into the LISA detection bandwidth. For quite 
generic initial conditions, such stellar clusters are expected to have reached a relaxed, 
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mass segregated, steady state which is independent of initial conditions at the time of 
formation. Strong (realistic) mass segregation is a robust outcome from the growth and 
evolution of stellar cusps around MBHs in the mass range 10 4 — 1O 7 M to which LISA 
will be sensitive. Our N-body results validate the Fokker-Planck description of the bulk 
properties of the stellar distribution. SMS boosts the EMRI event rates with respect 
to what would be implied by a shallower stellar density profile {e.g. 7 = 7/4, which 
has been the working assumption of almost all event rate estimates in the literature so 
far) that also respect the mass normalization obtained from observations of the Galactic 
center at 1 pc from the hole. In particular, our fiducial models of the Galactic center are 
enhanced by a factor of ~ 10 — leading to a predicted value of ~ 250 stellar bh EMRIs 
per Gyr. The FP formalism assumes two-body relaxation as the only dynamical driver 
present — this could be a severe restriction at radii < O.Olr^, inside which even the NB 
simulations with higher N in our sample start to run out o f particles and where RR could 
play an important role ( jHopman and Alexander! 120061 ; iMadigan et al.l |2010| ). Other 
crucial mechanisms are resonant relaxation, (small) triaxiali ty of the galactic potential , 



2007 



tidal separation of binaries and massive perturbers (see e.g. lAmaro-Seoane et al 
for a review). These are the subject of our current research work, and the extent to 
which they can significantly affect the EMRI rates is still an open question. 
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